Symmetry classification of spin-orbit coupled spinor Bose-Einstein condensates 
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We develop a symmetry classification scheme to find ground states of pseudo spin-1/2, spin-1, and 
spin-2 spin-orbit coupled spinor Bose-Einstein condensates, and show that as the S0(2) symmetry of 
simultaneous spin and space rotations is broken into discrete cyclic groups, various types of lattice 
structures emerge in the absence of a lattice potential, examples include two different kagaome 
lattices for pseudo spin-1/2 condensates and a nematic vortex lattice in which uniaxial and biaxial 
spin textures align alternatively for spin-2 condensates. For the pseudo spin-1/2 system, although 
mean-field states always break time-reversal symmetry, there exists a time-reversal invariant many- 
body ground state, which is fragmented and expected to be observed in a micro-condensate. 
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I. INTRODUCTION 

Ultracold quantum gases have provided an exception- 
ally idealized playground for emulating condensed matter 
systems |l|. Recently, Lin et al. have achieved a break- 
through in quantum simulations of condensed matter sys- 
tems by realizing an effective spin-orbit (SO) coupling in 
87 Rb atoms J2] . The SO coupling produces non-abclian 
gauge fields |3| , plays a key role in spintronics jj] such as 
spin-polarized transport, spin injection, and spin relax- 
ation, and yields many other interesting phenomena such 
as the quantum spin Hall effect and topological insula- 
tors [5j . The experimental realization of synthetic gauge 
fields has stimulated tremendous efforts on SO-couplcd 
quantum gases [6H18J. 

The SO interaction couples spin and linear momen- 
tum, thereby significantly modifying sin gle-p article spec- 
tra. For the Rashba-type SO coupling ^9l - |22j . the con- 
served quantity is L z + F z , where L z and F z are the pro- 
jected orbital and spin angular momenta along the z-axis, 
respectively. In homogeneous systems, an axisymmet- 
ric SO coupling will change the parabolic single-particle 
spectrum of a spin-F atom into 2F+ 1 energy bands with 
the lowest one featuring a Mexican hat, causing a circular 
degeneracy of single-particle ground states. For trapped 
systems, such circular dengeneracy is reduced into dou- 
ble and no degeneracy for half-integer and integer spin 
systems, respectively. 

In spinor Bose-Einstein condensates (BECs) [23|], the 
coupling between spin and linear momentum will cooper- 
ate or compete with spin-dependent or spin-independent 
interactions, giving rise to many exotic ground states 
with or without harmonic trapping potentials, such as 
plane-wave, stripe, triangular, and square lattice phases 
[a [ll|, UM ■ Firstly, the plane- wave phase and the stripe 
phase are found in SO-coupled pseudo spin-1/2 or spin-1 
BECs Q. Later, we found two different lattice phases 
where each spin component shows trianguar- or square- 
lattice density distributions in SO-couplcd spin-2 BECs 



with cyclic interactions [111 ]. As an axisymmetry har- 
monic trap is turned on, more phases are found as ground 
states of SO-coupled spin-1/2 BECs HIGH. Although 
many phases are found in SO-coupled spinor BECs, a 
systematic understanding on them is still elusive as pre- 
vious results largely depend on numerical solutions of the 
coupled Gross-Pitaevskii equations. In this article, al- 
ternatively we present a symmetry classification scheme 
to investigate ground states of SO-coupled spinor con- 
densates. We can then not only understand different 
lattice phases already found, but also find two differ- 
ent kagaome-lattice phases and a nematic vortex lattice 
phase, both of which emerge spontaneously without lat- 
tice potentials. 

This paper is organized as follows. Section II de- 
scribes the model Hamiltonian used in the present paper. 
Section III analyzes various phases that spontaneously 
emerge in SO-coupled spinor condensates based on sym- 
metry considerations. The following three sections dis- 
cuss properties of the three typical phases: triangular- 
, square-, and kagome-lattice phases. Section VII dis- 
cusses a possibility of a fragmented ground state of a 
SO-coupled spin-1/2 system with time-reversal symme- 
try. Section VIII summarizes the main results of this 
paper. 



II. MODEL HAMILTONIAN 

We consider a SO-coupled spinor BEC with TV atoms, 
including pseudo spin-1/2, spin-1, and spin-2 conden- 
sates in a pancake-shaped quasi-two-dimensional har- 
monic potential. The effective Hamiltonian is given by 
ti = Ho + 'Hint , with 



Ho — 

where ip = 
(x,y), V 
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(ipF, ■ ■ ■ , iP-f) t , M is the atomic mass, p = 
- Mu)\p 2 j2, v(> 0) describes the strength of 



the SO coupling, and F x>y are the spin-F matrices. 
For the F =1/2 case, we have 
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where g and g' denote the strengths of the intra- and 
inter-componcnt contact interactions, respectively. Here, 

«±i/2 = ^± 1/2 V'±i/2^±i/2V'±i/2, n = ni/ 2 + n_i/2, S z = 
n 1/2 - n_ 1/2 , a= (g + g')/2, and /3 = (g - g')/2. 

For the integer spin cases, we take the same interaction 
Hamiltonian from spinor BECs [23| ■ For the F = 1 case, 
the interaction part is given by 



H int , = ^ / dp[ cdltyfafa + MW k Pij ■ F k nPnP 3 ) , (3) 



where a and )3 give the strengths of density-density and 
spin-exchange interactions, respectively. Here, the in- 
dices that appear twice are to be summed over — F, . . . , F. 
For the F = 2 case, the interaction Hamiltonian is given 
by 
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Hmt = - / dp [aipl^^i + MWlFij ■ Fkrtrtj 
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where a, j3 and 7 give the strengths of the density- 
density, spin- exchange, and singlet-pairing interactions, 
respectively. Defining the harmonic-oscillator length 
a± = y/H/Mcj_i_, we introduce dimensionless parameters 
v' = v/u±a± and (a',/3',7') = (a,/3,'y)N/huj±. 



of the Hamiltonian, which in the present case is G = 
U(l) x SO(2) x T; (2) to list all subgroups H of G; 
(3) to find the order parameter that is invariant under 
H. If the order parameter can be uniquely determined 
from hip = tp for V/i € H, the state ifi is an inert state, 
which is always a stationary point of the energy func- 
tional. On the other hand, if the order parameter that 
is invariant under H is not uniquely determined, we need 
to minimize the energy functional within the restricted 
manifold. Such a state is called a non-inert state |25l - 
|28| . The high-symmetry states in the SO-coupled sys- 
tem are all non-inert states. To see this more clearly, 
we consider the eigenstates of simultaneous spin-space 
rotation C nzi where C nz = 1Zf(2it / 'n)lZp(2n / 'n) is the 
generator of a discrete cyclic subgroup of SO(2) with 
IZf = exp(— %F z 2-KJn) and 7?. p being respectively the 
2n/n spin and space rotation operators about the z- 
axis. For the integer-spin cases, C nz has n different eigen- 
values exp(— i2Nir/n) with N = 0, ...,n — 1, whereas 
for the F = 1/2 case, C nz has n different eigenvalues 
cxp{— i(2N+ l)ir/n} because (C nz ) n = — 1. If the eigen- 
value of C nz is not 1, we can infer that this state is invari- 
ant under the combined U(l) gauge transformation and 
spin-space rotation C nz , namely, O nz (N)ip = ip where 
CWN) = exp{i(2N + 1)tt /n}C nz for the F = 1/2 case, 
and O nz (N) = exp(i2Nn/n)C nz for the F = 1 and 2 
cases. For a spin-F 1 system, we can construct eigenstates 
of O nz (N) by using a complete basis set of plane waves 
as 
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III. SYMMETRY ANALYSIS 

Under a mean-field approximation, we assume that 
all atoms are condensed into a common single-particle 
state, and therefore the symmetries of the Hamiltonian 
are spontaneously broken. We can classify the grou nd 
states according to the remaining symmetries |24l428| . 

For spinor BECs [23( without SO-coupling, we usually 
focus only on the spin degrees of freedom because they 
are decoupled from the orbital degrees of freedom. As 
a result, the Hamiltonian is invariant under the U(l) 
global gauge transformation, the SO (2) spin rotation 
along z-axis for the F = 1/2 case or the SO(3) spin 
rotation for the integer-spin cases, and the time rever- 



sal T 



-iTT-F, 



H /C, where JC takes complex conjugation. 



In contrast, for SO-coupled spinor BECs, to make the 
Hamiltonian invariant, we should simultaneously rotate 
the spin and the space. Therefore, the Hamiltonian H 
described in Sec|H] is invariant under the global U(l) 
gauge transformation, simultaneous SO (2) global spin 
and space (spin-space) rotations, and time reversal T. 
According to the symmetry classification scheme [2j 
l28j . our task is: (1) to find the full symmetry group 



where k is the wave vector, S = {k| — tt/u + (p < tpf* < 
n/n + tp, (fk = arg(fca; + ik y )} with (p being arbitrary, 
\a) (a = —F, —F + 1, . . . , F) denotes the spin state with 
Alp = a under the z-axis quantization, and £>™(k) are 
the expansion coefficients. Unless we know the detail of 
£>™(k), we cannot uniquely determine the state that is 
invaraint under H = {E, O nz (N), . . . , OJJJ^N)} with E 
being the identity operator. Similar arguments can be 
applied to other subgroups of G. 

Since there are no inert states for spin-orbit coupled 
spinor BECs, we cannot find the order parameter by an- 
alyzing its symmetry only. However, such a difficulty can 
be alleviated if we take the ansatz that the state is a su- 
perposition of several degenerate single-particle ground 
states of the Hamiltonian in Eq. ((T|) with oj± = 0: 



V 



3=0 



VW(k = k g ,tp k = (p + j—), (6) 



where VW(k g ,ifk) = e lkcj ' p (-F(ipk) are degenerate 
single-particle ground states with k g = mv/h, ifk = 
aig(k x + iky), and Cp being arbitrary. For the cases of 
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FIG. 1: (Color online). Schematic diagrams for states preserving the combined symmetry of spin-space rotation C„ z , gauge 
transformation, and time reversal. Figures (a)-(f) correspond to n — 1, 2, 3, 4, 6, and oo, respectively, where arrows denote the 
direction of spin, solid closed loops indicate regions with nonzero momentum distributions |7/»(fc)| 2 , and n points on the circle 
are denoted as So, ■ ■ ■ , S n -i, respectively. Characteristic momentum distributions of the ground state ip for a trapped system 
are shown in the lower right corner of each figure. Black (yellow) color refers to the region of small (large) amplitude. 



F = 1/2, 1 and 2, we have respectively 

C-i/2(^) = (l,-e^) T /V2, 
C-i(^) = (e-^ fc ,-\/2,e^) T A 
C- 2 Ofc) = (e- 2 ^*,-2e- i ^, V6,-2e^,e 2l ^) T /4-(7) 

The corresponding spin expectation value is anti-parallel 
to the wave vector: ((F x ),(F y )) = — F(cosf k, sin ip^). 
If we take <j>j = j2Ntv/ti, the state ijj in Eq. (JB|) is in- 
variant under O nz (N). This state is a special case of 
that in Eq. (J5} with the set S shrinks into only one el- 
ement and the allowed values of k are fixed at discrete 
points as -D™(k) = S(k — k g )[(-F(<p)]o-- Figure 1 illus- 
trates schematic diagrams for the states of Eq. (JSJ) , where 
n single-particle ground states PW(k g ,fk) are denoted 
by points Sj. 

We choose the ansatz in Eq. ([5]) based on previous 
understanding on a 2D system and our numerical results 
given in the following sections. For a 2D homogenous sys- 
tem with uj_ = 0, as discussed in Rcfs. [s[ll|, plane- wave 
and stripe phases are found to be the ground states for 
the F = 1/2, 1 and 2 condensates @, HH, whereas trian- 
gular and square lattice phases exist only for spin-2 con- 
densates with cyclic interactions [11] . Such phases can 
all be described by Eq. © . For a trapped system, trans- 
lation symmetry is broken, and linear momentum is no 
longer conserved. Therefore, in momentum space, there 
arc infinite points forming regions, where momentum dis- 
tributions |-0(k)| 2 of the ground state ip are nonzero. In 
Fig. [TJ we use solid closed loops (circles for n — oo ) to 



denote such regions. In the case of strong SO couplings, 
we can approximate low-lying single-particle eigenstates 
as HHl 



*,- 



v {k') oc fc'-V 2 e -(*'-*') 2 /2#_.( fe ' 



v')e im * k (MVk),(8) 



with eigenenergies En >m = (m + l/2) 2 /2t/ 



72, (m 2 + l/4)/2u' 2 + n + 1/2 - 
2 /2 for F = 



n + 1/2 - 
72, and (m 2 + 
3/4)/2i/ 2 + n + 1/2 - v' 2 /2 for F =-- 1/2, 1 and 2 re- 
spectively. Here, n = 0,1,2,..., m = 0, ±1,±2, . . ., 
k' = ka±, and H^ are Hermite polynomials. When the 
interatomic interactions are weak, the ground states can 
be constructed only from the states with n = 0, which 
are superpositions of single-particle eigenstates £-F(k) 
with weight j fc /-V2 e -(fc'-f')V2 ~ v /-i/2 e -(fc'-«')V2. In 
the limit of w_l -¥ 0, e -(fc'-»') 2 /2 _> S(k'-v')/V2lr, which 
implies that the ground states will involve momenta with 
its magnitude close to k g . Furthermore, by choosing spe- 
cific weak interactions, axial rotational symmetry in the 
single-particle states of Eq. ((8]) can break into discrete 
rotation symmetry, and 2D latticephases appear, as al- 
ready found in the F = 1/2 case [lj, [Taj. In the limit of 
weak trapping potential, strong SO couplings, and weak 
interatomic interactions, there are ground states whose 
momentum distributions split into several small regions 
as those illustrated in Fig. [TJ As long as the area of the 
allowed momentum is small enough, these ground states 
can still be well described by Eq. ([6]). 

The validity of the ansatz has been numerically con- 
firmed [30j, because we find several ground states which 
can be understood by Eq. ^. In the lower right corner 
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FIG. 2: (Color online). Ground states of trapped SO-coupled spin-1 condensates, with v = 15, a = 0.5 and /?' = 0.1. (a) 
Density distributions of three spin components Mf — 1,0,-1 from left to right, (b) Spatial variation of the corresponding 



order parameter visualized by plotting £^ 
of the spin wave function. 



£ m Y F n (0, ifi), where Y F F is the spherical harmonics and £ m is the m-th component 



of each figure in Fig. [TJ we illustrate the characteris- 
tic momentum distributions of these ground states. For 
the case of n = 1 and 2, the state in Eq. ((SJ) can de- 
scribe respectively the plane-wave phase and the stripe 
phase found in a homogeneous system [8|, [TlJ . For the 
case of n = 00, the state is invariant under the combined 
SO(2) spin-space rotations and U(l) gauge transforma- 
tion, and has been predicted for the F = 1/2 case [bffla ]. 
where single-particle eigenstates of Eq. (|5J) will be stabi- 
lized in some parameter regions. For the F = 1 and 
F = 2 cases, similar arguments can be applied, and more 
single-particle eigenstates of Eq. ([5J can be stabilized. 

In the following, we start from the ansatz of Eq. ^ . 
and find several high-symmetry states, where a set of fa 's 
are determined by requiring the wave function of Eq. © 
to preserve a certain symmetry. Some of them become 
ground states of a trapped SO-coupled system, where 
triangular, square and kagaome lattices are formed spon- 
taneously without lattice potentials. These states can be 
described by Eq. ([6]) with n = 3, 4, and 6, respectively. 

Before discussing lattice phases, we point out that (1) 
the spin and space rotation operations commute with 
time reversal [T,C nz ] = 0, and that (2) for the F = 1/2 
case, the ground-state order parameter ip always breaks 
time-reversal symmetry because T 2 = — 1, which means 
ip will only be invariant under C nz or TC nz , and possibly 
combined with the U(l) gauge transformation for each 
n. For the integer-spin case, the ground-state order pa- 



rameter will be invariant under one or all of the three 
operations: T, C nz and TC nz , and possibly combined 
with the U(l) gauge transformation for each n. 



IV. TRIANGULAR-LATTICE PHASE 

A triangular-lattice phase is described by Eq. © with 
11 = 3 as 

V-> = e^°e lfcs:E C-F(0) + e l * 1 e^ ( - a;/2+v ^ y / 2 )C-F(27r/3) 
+e ^e lk ^ x/2 -^ y/2) t- F (An/?>), (9) 

where we take (p = as an example. According to 
the symmetry classification scheme described in Sec lIIH 
we can apply specific symmetries on the state of Eq. 
© to evaluate the value of fa. We find the only 
possible symmetry is described by the group H = 
{E,0 3z (N),0| z (N)}, resulting in fa = j2Ntt/3 + const. 
Meanwhile, the states of Eq. ((9j) with different values 
of fa are connected by a global U(l) phase change and 
a global coordinate translation. This is because for ar- 
bitrary change 5 fa in fa, there are always solutions for 
(6x, 5y) and 8U that satisfy 

cxp{i[5(po + k g Sx + 5U)} = 1, 
exp{i[8fa + k g (-Sx/2 + VS5y/2) + 5U)} = 1, 
cxp{i[(50 2 + k g (-Sx/2 - VZSy/2) + 5U)} = 1, (10) 



where (5x, Sy) and SU describe the amount of the coordi- 
nate translation and the U(l) phase change, respectively. 
As a result, there is only one type of triangular-lattice 
phase. By choosing a proper center of lattice, such a 
ground state becomes invariant under 3z . Note that in 
a harmonic trap, the 3z symmetry axis sometimes devi- 
ates from the trap center to gain the interaction energies. 

In Fig. [2] we illustrate a triangular-lattice phase which 
appears in the ground state of SO coupled spin-1 con- 



shift, the state i/> in Eq. (JTTJ) changes with 



densates with 



15, a' = 0.5 and j3' = 0.1, showing 



triangular-lattice density distributions for each spin com- 
ponent. Figure [2jb) shows the spatial variation of the 
corresponding order parameter. At three lattice sites, 
the order paramters can be written as V\ : (0, e~ l7r / 3 , 0) T , 
V 2 : (0,e i7r / 3 ,0) T and T 3 : (0, e l7r ,0) T . To go from Vi to 
V2, a 7r spin rotation along V\Pi and the e~ I7r / 3 gauge 
transformation are needed. Similar transformations are 
needed to go from V2 to V$ or from V3 to V\ . Going along 
the loop oiV\P2VzP\ 1 the order parameter undergoes a 
■k spin rotation along an axis on the x-y plane and the 
— 7r gauge transformation, which, however, does not im- 
ply that the mass circulation here is fractional with —1/2 
winding [31(. To check whether this vortex is fractional 
or not, we should calculate the circulation [23[ along the 
loop of V-lP-xVzPx. Numerically we found this vortex is 
not fractional. 

For the triangular-lattice phase, two (three) spin com- 
ponents of spin-1/2 (spin-1) condensates tend to occupy 
different spaces due to the interferences of three plane 
waves. For spin- 2 condensates, five spin components 
are divided into three classes: (1) Mp = 2,-1; (2) 
Mp = 1, -2; (3) Mp = 0. In the cases of (1) and (2), two 
different spin components show the same wave function, 
which can be understood from Eq. ([9]) . These properties 
are consistent with our previous work |ll| . 



V. SQUARE-LATTICE PHASE 

When n — 4, the state in Eq. ^ can be simplified as 

V> = e # V fc s x C-F(0)+e 4 *V fc ^C-F(7T/2) 

+ e*e-*«*(.F()r) + e^ 3 e~ tfc ^C-F(37r/2), (11) 

where we again take <p = as an example. This state can 
only potentially be invariant under O nz with n = 1,2,4. 
If the state is only invariant under 0\ z or 2z , V' cannot 
be uniquely determined. Therefore, we require the state 
preserving 04 2 (N), resulting in fa = jNir/2 + const. 

Similar to the triangular-lattice phase, there are other 
states which are related to the states invariant under 
04z(N) by simply doing a global lattice shift. They may 
be classified into the same class, as they have the same 
lattice structure. The criterion for each class is deter- 
mined by the parameter A = (fa + fa) — (fa + fa). The 
reason is given as follow: When doing a global lattice 
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(12) 



Absorbing Sx and Sy into fa , we find that there are two 
invariants fa + fa and fa+fa. Their difference A is also 
invariant under the global U(l) gauge transformation. 

Furthermore, we find that for the states invariant un- 
der 4z (N), A = (2M + 1)tt when N = 1,3; A = 2Mtt 
when N = 0, 2, (M g Z). We then define the criterion 
for two different types of square-lattice phases with dif- 
ferent lattice structures as (1) A = (2M + l)ir and (2) 
A = 2M7T. The states that fall into the same class are 
connected by a global lattice shift. By doing a proper lat- 
tice shift, the symmetry of the state tp can be described 
by H = {E 7 4z (N),0 2 4z (N),0 3 4z (N)} with N = 1,3 for 
the first class and N = 0, 2 for the second class. For 
the integer-spin cases, the state ip in the second class 
is also invariant under time reversal T |32| by choosing 
a proper lattice center. Numerically, we find all such 
ground states to be integer-spin condensates. Again, in 
a harmonic trap, the 4z symmetry axis does not always 
coincide with the trap center to minimize the interaction 
energy. 

The ground states for the case of A = (2M + l)7r are 
found to exist in spin-1 BECs with j3 < 0, and spin-2 
BECs with (3 > and 7 > jjl|. Figures Ufa) and 
(b) show the corresponding ground-state density distri- 
butions and order parameters for the spin-1 BECs, with 
v' = 15, a' = 0.5, and /?' = —0.1. The spin compo- 
nent Mp — forms a square lattice with lattice con- 
stant 7r//c g , while the spin components Mp = 1 and —1 
fill the center of squares alternatively, both forming a 
square lattice with lattice constant \/2tt /k g . The numer- 
ically obtained total density distributions are smooth for 
such ground states, consistent with the prediction from 
the corresponding state tp m Eq. ([6]). Only the compo- 
nent with Mp = involves vortices, with vortex cores 
filled by the Mp = 1 and Mp = — 1 components having 
vorticity H and —h, respectively. Thus, the spin vortex 
lattice is filled by ferromagnetic vortex cores which are 
polarized in the z and —z directions alternatively. For 
the F = 2 case, such phase has already been predicted in 
our previous work [ll| . We have found that spin compo- 
nents Mp = 2, 0, —2 show the same density distribution 
which is different from those of the Mp — 1 and —1 com- 
ponents. Again, there is no fractional or integer vortex 

Another type of square lattice phase with A = 2M-7T 
appears in the trapped SO-coupled spin-2 condensates 
with antiferromagnetic spin-dependent interactions. Fig- 
ure 02c) shows ground-state density distributions of spin 
components Mp = 2, 1, from left to right, with v' = 15, 
a' = 0.5, j3' = 0.1, and 7' = —0.1. The spin components 
Mp = — 1 and —2 show the same density distributions as 
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FIG. 3: (Color online). Ground states of trapped SO-coupled spinor BECs, with v — 15. (a) Density distributions of Mf — 1 
(left), (middle), —1 (right) of the spin-1 case with a' = 0.5, /3' = —0.1. (b) The corresponding order parameter, where F is 
used to denote ferromagnetic vortex cores, (c) Density distributions of Mf = 2 (left), 1 (middel), and (right) of the spin-2 
case with a' — 0.5, j3' = 0.1 and 7' = —0.1. (d) The corresponding order parameter. The uniaxial and biaxial nematic vortex 
cores are denoted by UN and BN, respectivley. Here, black (yellow) color refers to the low (high) density region. 



Mf = 1 and 2, respectively. In this phase, we find that 
two physical quantities \{F)\ = \Y,ij ^Fij^A/ T^i^i^i 
and |(e)| = I Eiji-iy^^r-i^-jl/Ei^^i areuni- 
form and are equal to and 1, respectively. Figure EJd) 



illustrates the corresponding order parameter, showing a 
spin vortex lattice with uniaxial nematic (UN) and biax- 
ial nematic (BN) vortex cores aligning alternatively. 




FIG. 4: (Color online). Ground states of trapped SO-coupled pseudo spin-1/2 BECs, with v' = 15 and (a,b) a' — 0.3, j3' — 0.06, 
(c,d) a' = 0.4, f3' = 0.08. Here, black (yellow) color refers to the low (high) density region. From left to right, we show the 
spin-up, spin-down and total density distributions. Figures (a) and (c) are numerically calculated for a trapped system, while 
(b) and (d) are obtained by superposition of six single-particle states VW(ipk)- In (b) and (d), we use symbols "-","+", and 
"++" to denote vorticities —h, h and 2ft, respectively. The size of figure (b) [(c)] is the same as that of (d) [(a)]. 



VI. KAGOME-LATTICE PHASE 

The ansatz in Eq. (J5J) with n = 6 can be written as 

= e #o e lfc « x C-F(0)+e #1 e 4 ^ (;E/2+x/5!//2) C-F(7r/3) 
+ e i ^e ik ^- x/2+VSv/2) C- F (2n/3) + e ilp3 e- %k * x C~F{K) 
+ e l * 4 e ife9( " ;r/2 "^ y/2) C-F(47r/3) 

+ e ^ e ^ S (^/2-V / 3y/2) C _ jF ( 57r / 3 ) j ( 13 ) 

where we take tp = as an example. This state can 
describe kagaome-lattice structures. Similar to the case 
of n = 4, we note that four independent phases 

A = 00 + 02 + 04, Ai = 00 + 03, 

A 2 =02+05, A 3 =0 4 +01, (14) 



are invariant under a global shift of the lattice. Besides 
the one used to describe the global phase change, there 
will be three independent phases which determine the 
structure of the lattice. Therefore, there are infinite types 
of lattice structures of the state in Eq. (fl~3")) . 

According to symmetry classification scheme, we need 
to apply specific symmetries to determine the state ip. 
Only if the symmetry is high enough, the state ip can be 
uniquely determined. For the F = 1/2 case, two high 
symmetry groups are generated by operators Oq z and 
TOq z , respectively, whereas for the integer-spin cases, 
high symmetry groups can be generated by only Oq z or 
{0 6z: T}. 

(1) If we require the state is invariant under 06z(N), 
the state ip is uniquely determined with <pj = jNtt/3 + 
const. For the integer-spin cases, the state is further 



invariant under time reversal T [32| if N = 0,3. 

(2) For the F = 1/2 case, there are states that pre- 
serve the TOq z symmetry and break the 062 symmetry. 
To determine such states, we start from the states in- 
variant under C 3z (N) with N = 0,1,2. The value of 
<f>j is determined as <j>j = j2N7r/6 + mod(j, 2)<f> + const, 
where </> is arbitrary. Furthermore, by requiring that the 
state be invariant under TC^z |32| , we obtain N = 1 and 
(j) = tt/2 + Mtt. 

We numerically find two distinct classes of kagaome- 
lattice phases in a spin-1/2 BEC with j3' > 0, whose 
symmetries are described by the group H generated re- 
spectively by (1) cxp(i27r/3)C3 Z and TCi z (with <pj = 
J7r/3+mod(j, 2)(l/2 + M)tt + const) and (2) exp{i(M + 
l/3)ir}C ez (with fa = j(M + 1/3)tt + const), where the 
ground states with different values of exp(iM7r) in the 
same class are time reversal with each other. Figure Q] 
shows the spin-up, spin-down and total density distri- 
butions for numerically obtained order parameters in a 
trapped system [Fig. d£a) and (c)] and the correspond- 
ing ansatz of Eq. (fT3"]) [Fig. |3Jb) and (d)]. For each 
spin component, density distributions show kagomc lat- 
tice structures. 



VII. FRAGMENTED GROUND STATES 

Although we have discussed within the mean-field the- 
ory so far, there is a possibility that a fragmented ground 
state [33H3q |. rather than a mean-field state, arises. Ac- 
tually, for the case of SO-coupled spin-1/2 system with an 
SU(2) symmetry, we construct a fragmented state whose 
energy is degenerate with the mean-field solution up to 
the mean-field approximation. Such a fragmented ground 
state is expected to arise in a mesoscopic system. In this 
section, we discuss fragmented ground states from the 
point of view of the time-reversal symmetry. 

As pointed out in the Scc lIIIl mean-field states always 
break time reversal for the F = 1/2 case because T 2 = 
— 1. In this section, we show that there are also time- 
reversal invariant many-body states that are fragmented 
in a time-reversal preserving F — 1/2 bosonic system 
with or without SOCs. 

For a pseudo spin-1/2 system with N atoms, we have 
T 2 = (-) Ar . When N is odd, all states break time- 
reversal symmetry. In contrast, when TV is even, there 
is always a time-reversal invariant ground state which is 
fragmented: if [T,H] = and H\ip) = E\ip), \ip T ) = 
\ip) + T\ip) is an cigenstate of H with the same eigenen- 
crgy E, and is invariant under time reversal if N is even. 
To check whether \ipj-) is fragmented or not, we can di- 
agonalize its single-particle density matrix p [36J. Wc 
can regroup the single-particle eigenstates as { v I / Mi , ^ Vi }, 
where we use \x and v to distinguish two states by their 

-* Mi . Time- 



the states { x I / Mi , ^ Vi } are 



time reversal as T^ • 



^ Vi and T*„< 



reversal properties of the corresponding creation opera- 
tors {&'„ ,al} and annihilation operators {0^,0,^} for 



Ta u r- 1 



T 



TalT- 1 : 



(15) 



We also have {ipT\a} li a tlj \ip T ) = (TV'rlTa^a^ |Vr)* = 
(tjjj-\al,.a v - \iI>t)* i where for the first equality we 
use the fact that T is an antiunitary operator, 
while the second equality is due to Eq. I|15p 

and the fact that \ipf) is invariant under time re- 
versal. Similarly, we obtain (ip'j-\a'l.a I y.\ipj-) = 

(^rlal^lipr)* > (V'rlo^a^lVr) = -(^rja^a^l^r)*, 
and (i/j T \al.a H \ilj T ) = -(ip T \aj lz a^\i> r )* ■ These equal- 
ities imply that the single-particle density matrix is in- 
variant under S as 



■p~-i 



Ss jcn< 



1 
-1 



(16) 



Due to this symmetry, we can infer that if |\&) is an cigen- 
state of p, S|W) is also an eigenstate of p with the same 
cigenencrgy but orthorgonal to |W). Our argument above 
is also valid for the ground state. 

A simple state that preserves time-reversal symmetry 
is 



1 / \ C-t \ N/2 (-1 \ N/2 



vac) 



(17) 



Using this ansatz wave function, we numerically mini- 
mize the energy functional {ipq-\T-L\ipq-) , and find that the 
ground-state energy is the same as that obtained with 
single-orbital mean-field approximation for an SU(2)- 
symmetric system without spin-dependent interactions, 
where we take the lowest-energy band approximation by 
considering only states \&o,z with v' = 15, f3' = 0, and a' 
ranging from to 1. 

Before conclusions, we would like to point out that 
the fragmented ground states we discuss here are pro- 
tected by time-reversal symemtry for a spin-1/2 system. 
The two-orbital ground states for an SU(2)-symmctric 
spin-1/2 system we constructed in Eq. (fT7|) arc different 
from that predicted in Refs. [33, [38| in which the frag- 
mented ground states are produced in a cooling process 
that conserves the total spin S. For an initial state at 
high temperature, we have S ~ vTV. If all atoms con- 
dense into the same orbital state, S = N/2. Therefore, 
to be compatible with the requirement of the total spin 
conservation, atoms should condense at least into two 
orbitals, resulting in a fragmented ground state. 



VIII. SUMMARY 

We have systematically classified ground states of 
strong SO-coupled trapped spinor BECs, including 
pseudo spin-1/2, spin-1 and spin-2 cases, based on sym- 
metry analysis. In accordance with breaking of simul- 
taneous SO(2) spin-space rotation symmetry in favor of 



discrete symmetries, there emerge lattice phases showing 
stripe, triangular, square, and kagome lattice structures 
on each spin component. Imposing symmetries of T, C nz 
or TC nz or combined them with the U(l) gauge transfor- 
mation on the order parameter, we predict several lattice 
phases, some of which are found to be ground states of 
a trapped system. For the spin-1/2 case, the ground 
states can be classified into two classes: one breaks time- 
reversal symmetry and the other is invariant under time 
reversal and shows fragmented properties. 
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